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ABSTRACT 

The study of asteroid families has provided tremendous insight into the forces that sculpted the main belt and 
continue to drive the collisional and dynamical evolution of asteroids. The identification of asteroid families 
within the NEO population could provide a similar boon to s tudies of their form ation and interiors. In this study 
we examine the purported identification of NEO families bv lDrummondl J200(!) and conclude that it is unlikely 
that they are anything more than random fluctuations in the distribution of NEO osculating orbital elements. 
We arrive at this conclusion after examining the expected formation rate of NEO families, the identification 
of NEO groups in synthetic po pulations that conta in no genetically related NEOs, the orbital evolution of the 
largest association identified bv lDrummondl J2000I) . and the decoherence of synthetic NEO families intended to 
reproduce the observed members of the same association. These studies allowed us to identify a new criterion 
that can be used to select real NEO families for further study in future analyses, based on the ratio of the number 
of pairs and the size of strings to the number of objects in an identified association. 
Subject headings: asteroids, dynamics 



1. INTRODUCTION 

The discovery of a group of Near Earth Objects (NEO) 
that all derive from the products of a catastrophically dis- 
rupted NEO parent body would be an important step in un- 
derstanding their structure and dynamical evolution. It would 
be equivalent to the information afforded a geologist in crack- 
ing open a rock with a hammer. Studies of the chips off the 
old block of a larger asteroid will provide interesting insights, 
just as the study of similarly formed main belt (MB) aster- 
oid families has pro vided critical information to the study o f 
those asteroids (e.g.. lCellino et al.l2004HZappala et al.l2002l) . 
The difference is that NEO families will have a much faster 
dynamical evolution such that they are detectable over much 
shorter timescales as they diffuse into the orbital element 
space of the background objects. They also sample a much 
different range of parent asteroid sizes since the known sam- 
ple of NEOs are much smaller than the known sample of MB 
objects. 

The time scale of catastrophic collisions between MB as- 
teroids that create asteroid families is short compared to the 
age of the Solar System. In catastrophic collisions, the col- 
liding bodies completely shatter into smaller fragments, some 
of which are dispersed and then re-accumulate into objects on 
independent but similar orbits dDavis et alJl2002l) . The first 
groups of aster oids with nearly th e same orbital elements were 
discovered by iHiravamal ( 119181) and these groups are now 
known as Hirayama families. The members of each family 
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usually share similar spectral properties dCellino et al.l2002l) . 
further suggesting a common origin in a single parent body 
that has undergone a catastrophic collision. We consider the 
members of this kind of asteroid family "genetically" related 
since they share the same parent body. At this time, more than 
50 significant familie s have been tabulated in the MB (e.g., 
iNesvornv et al.l2005l) . 

MB families are identified by searching for enhancements 
in the space of proper orbital elemen ts (a„, e„, A prec ise 
definition of the proper elements (e.g. lKnezevic et al.l2002l) is 
beyond the scope of this paper but it is sufficient here to think 
of them as long-term average orbital elements. This method 
is not readily applicable to NEOs for two main reasons: (1) 
the calculation of proper elements for NEOs is problematic 
due to the extremely chaotic evolution of NEO orbits and, (2) 
the limited population (3,319 objects as of 2005 Apr. 28) of 
known NEOs are spread over a much larger orbital element 
space than MB asteroids. Therefore, while asteroid families 
have long been recognized in the MB, the situation remains 
unclear regarding the existence of genetically related families 
other than meteor streams in the NEO population. 

We have identified five possible mechanisms for the pro- 
duction of genetic NEO families: 

1 . NEO-MB asteroid collisions, 

2. tidal disruption of NEOs, 

3. intra-NEO collisions, 

4. a MB family producing event near a MB resonance 
followed by rapid dynamical evolution of some of the 
members into NEO orbits, and 
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5. spontaneous NEO disruption. 

Among the five, we consider the first two to be the most 
likely method for producing NEO families. Collision s be- 
tween NEOs are expected to be rare ( iBottke et al!l996l) since 
the density of NEOs interior to the MB is miniscule compared 
to the density of asteroids in the MB. The fourth mechanism 
is unlikely because there is only a small fraction of MB or- 
bital element space near resonances capable of transporting 
the residue onto NEO orbits. Furthermore, the chaotic dynam- 
ical evolution during the transportation of the fragments to 
near-Earth orbits makes it unlikely that such a family could be 
identified as such based on their osculating elements. The last 
mechanism is reserved for objects with high volatile content, 
presumably comets that have worked their way i nto NEO or- 
bits, a nd is probably the way meteor streams (e. s. . iDrummondl 
are created. Thus, meteor streams are genetic NEO 
families but the members are orders of magnitude smaller than 
the smallest detectable NEOs. 

The number of catastrophic NEO collisions per year that 
may produce families due to collisions on targets up to abso- 
lute magnitude H max is given by: 

/H m a x 
n(H)p c (H)dH (1) 
-oo 

where n(H) is the differential number distribution of 
NEOs, and pc{H) is the annual probability of catastrophic 
collision, both given as a function of absolute magni- 
tude. Using the NE O size-frequency distribution (SFD) of 
IBottke et ail <2002al) and the collision probabili ty for MB 
asteroids from Fig. 14 of IBottke et al.l (120051) . we esti- 
mate that during the typical 10 6 year dynamical life time 
of NEOs dMorbidellietaUl2003 iGladman et all 119971) <?SC 
l/~8/~1000 NEOs larger than 10/1/0.1 km diameter are 
catastrophically disrupted. However, the catastrophic col- 
lision probability for NEOs is dramatically smaller than 
that of MB asteroid s due to their much lower space-density 
jBottkeetalJlT996l) . Once an NEO is dynamically decou- 
pled from the MB the collision probability drops even further. 
Thus, we consider it unlikely that any NEOs > 1 km diameter 
have suffered a family-producing event while the number of 
NEOs >100 m diameter that have produced families is small. 

An asteroid family produced in the disruption of a target 
would produce a large number of smaller asteroids. The SFD 
of the re-accumulated frag ments may be determined through 
nume rical modelling (e.g., iMichel et alJ 120021 iDurda etafl 
120051) or through studies of the SF D of known MB aster- 
oid families (e.g.. lTanga et al.lfl999l) . The SFDs of the fam- 
ilies are usually described by a power-law in the diameter 
D, N(> D) oc D-p, where 2 < p < 5 for real MB 
families (e.g., Eunomia , Flora & Koronis) though it appears 
JMorbidelli et all2003l) that the SFD of observed families be- 
comes shallower for H > 15 (D < 5 km). The SFD of 
simulated new asteroid families displa ys a rich morphology 
JMichel et alJl2002t IDurda et ai1l2005l) with a wide range of 
slopes that is difficult to characterize with a single power-law 
over the entire range of sizes but they are qualitatively similar 
to the observed SFD. 

For the purpose of quickly estimating the expected number 
of known fragments due to a NEO family producing event, we 
assume that the fragments have a SFD similar to the known 
families with 2 < p < 5. In a barely catastrophic collision the 
largest fragment (LF) has half the mass of the target. Since 
we know that M oc V oc D 3 and that the diameter of an 



asteroid is related to its absolute magnitude by D oc 1Q~ H ' 5 , 
then it is simple to show that Hlf = Ht + 0.5 where Ht 
is the absolute magnitude of the target. By definition, N(< 
H LF ) = 1, and thus N(< H) = \q^h-h t -o.5) so that me 
differential HFD for the fragments must be: 

n(H) ~ 2.3 a iq^h~h t -o.5) (2) 

with 0.4 < a < 1.0 and Ht being the absolute magnitude of 
the target. 

If the actual and known population in a volume dadedidH 
are represented by dN and dn respectively, then dn = 
e(a,e,i, H)dN where e is the observational completion of 
NEOs with orbital elements in the range a — > a + da, 
e — > e + de, i — > i + di and H — > H + dH. The number of 
known fragments from the collision is then 

/+oo 
n(H) e(H) dH (3) 

where e(H ) ~ e(a', e', i', H) and (a', e', i') are the mean or- 
bital elements corresponding to a particular family. 

Using reasonable estimates for e(H) and the full range of 
possible slopes for the SFD, we find that the number of de- 
tectable fragments in a NEO family produced in the collision 
with a 1 km diameter target varies from essentially zero to 
many tens of thousands (if the family has the most steep SFD). 
This has tremendous implications for the formation of haz- 
ardous NEO families and spikes in the collision probability 
with the Earth. However, as argued above, the catastrophic 
disruption of a 1 km diameter NEO is highly unlikely. In 
the case of a 100m diameter NEO target the number of ob- 
servable fragments is essentially zero. The fragments would 
be very small and would quickly begin to evolve apart dy- 
namically and be strongly a ffected by non-gravi tational forces 
(e.g., Yarkovsky forces: see lBottke et all2002bL for a review). 
Our ability to recognize a genetically related NEO family will 
depend on the efficiency of the employed algorithm and the 
time scale for decoherence of the members" orbits. We will 
explore both these factors in this work. 

Thus, it would seem unlikely that observable NEO fami- 
lies can be produced and identified in the known NEO pop- 
ulation. That being said, the opportunity of finding an NEO 
family would be important to asteroid physical studies and 
would also allow us to refine our understanding of the catas- 
trophic disruption p r ocess and statistics. In a similar vein, 
iPauls and Gladmanl ll2005t) have studied the decoherence of 
"meteorite streams" to determine if time-correlated meteorite 
fa lls exist. 

iDrummo'ndl d!99 ll 120001) made the first attempts at search- 
ing for groupings in the known NEO population. His study 
was based on a £>-criterion analysis (orbit similarity) of oscu- 
lating orbital elements. Osculating elements are the instanta- 
neous 2-body (Sun+object) orbital elements for an object. He 
discovered 14 associations, 8 strings and 7 pairs out of a sam- 
ple of 708 NEOs. Briefly, associations are density enhance- 
ments above the local background in a four-dimensional D- 
criterion space (the values in each of the four dimensions are 
functions of an object's five osculating elements not includ- 
ing the mean anomaly). Strings are sets of NEOs in which ev- 
ery member is "connected" to at least one other member with 
D-criterion less than a specified threshold value, and pairs are 
two objects that have orbits that are statistically determined to 
be improbably close. These groupings contain a total of 155 
NEOs or ~ 22% of the entire sample known at the time of 
Drummond's study. 
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At this time, 56 objects in his groupings (36%) have been 
taxonomically classified 1 but there is no obvious correlation 
between the taxonomic types in his groupings. Moreover, 
there are cases where very different taxonomic types exist 
within the same grouping — for instance, both S-type and 
F-type NEOs exist in the Al a s sociat ion (the largest associ- 
ation identified by iDrummondl (I2000I) . see Tabled). In MB 
families there is a high degree of correlation between the col- 
ors o r taxonomic types of objects associated wi th each family 
(e.g- lNesvornv et alJl2005llCellino et al.l2002l) . The fact that 
there exists heterogeneity in Drummond's (2000) NEO asso- 
ciations may be used to argue that the groups are not geneti- 
cally related families. Alternatively, it may be argued that the 
groups represent the fragments of an inhomogeneous parent 
body. The latter explanation is probably unlikely given that 
the NEOs are themselves probably fragments from collisions 
that took place within the MB. 

In this paper we test the significance of the purported 
NEO families through various numerical simulations. First, 
we briefly describe (jffl the analysis technique employed by 
lDrummond ( n^9lll2000i) . In we discuss the statistical sig- 
nificance of the detections by attempting to search for group- 
ings in synthetic NEO populations where no genetic families 
exist. We then demonstrate the decoherence of synthetic Al- 
like families as they evolve under only the effect of gravita- 
tional forces in 0and study the dynamical evolution of the 
actual Al association in §5\ These studies allowed us to de- 
velop a technique for identifying real families within the NEO 
population as described in i|6] We discuss our main results in 
SQand close with a summary in ||8] 

2. NEO FAMILY IDENTIFICATION 

Several ve rsions of the D-criterion have been intro- 
duced (e.g., ISouthworth and Hawkins I Il963t IDrummondl 
H98ll Ijopekl Il993t IValsecchi et al] 119991) to quantify or- 
bital simila r ities. The dimensionless D-criterion used in 
|Drummond| (| 2000|) and adopted in this study was defined by 
ISouthworth and Hawkins I ( 119631) as: 

D = ^df + 4+dl + dl (4) 

where 

72 ill ~ <?2 ^2 ,2 / \2 

1 = ( AU ' ' 2 = ^ ~ E2 ' ' 
d 2 = [2sm{I/2)] 2 ; d\ = [(d + e 2 ) sin(n/2)] 2 (5) 

and 

/ = arccosfcosi! cosi 2 + sin^ sini 2 cos^! — f2 2 )] (6) 

TT i o • r U + «2 . ^i - ^2 !i 

11 = lo\ — UJ2 + 2, arcsm[cos — - — sin sec —J. (7) 

The subscripts 1 and 2 refer to the two orbits being compared, 
and the sign of the arcsin term of Eq. ^changes when | fl\ — 
il 2 |> 180°. q = a(l — e) is the perihelion distance, e the 
eccentricity, i the inclination, co the argument of perihelion, 
and £1 the longitude of the ascending node. 

The £>-criterion is a dimensionless distance metric in a 
four-dim ensional space. Th e volume search method intro- 
duced bv lDrummondl ( 12000 ) to identify associations searches 
for density enhancements above the local background around 
every object in the sample in this four-dimensional space. It 

Data collected from the PDS Asteroid/Dust Subnode, 
http://www.psi.edu/astsubnode.html 



estimates the local background level by fitting a uniform lo- 
cal background and a Gaussian association (a gamma distri- 
bution) simultaneously to the cumulative distribution of D 2 
around each orbit 2 , and a candidate group is identified if the 
association plus background model gives a better fit to the data 
than the background-only model. Those asteroids with D- 
criterion less than the value at which the association density 
equals the background density are gathered to form a mean 
orbit around which the cumulative distribution of D 2 is fit for 
the two models again. The procedure iterates until it con- 
verges on a stable group of objects. Th e groupings detect ed 
in this way were termed associations b v IDrummondl d2000l) . 

We also searched for pairs and strings base d on the D- 
criterion in a manner similar to that employed bv lDrummondl 
J2000l A couple of NEOs are a pair if the £>-criterion be- 
tween the two is below a cutoff D pa i r . Similarly, a string of 
NEOs is a set of objects in which every member is connected 
to at l east one other with a D < D str ing- In IDrummondl 
(2000), a cutoff of D st ring = 0.115 was chosen under the 
guidance of the results from the volume search. We set a more 
stringent threshold, D str ing = D pa i r — 0.1, throughout this 
paper for the reason discussed below (fj6ji. 

In Drummond's (2000) study, the search for groupings was 
performed on a sample of 708 NEOs (through 1999 HF1) 
listed on the Minor Planet Center (MPC) homepage as of 1999 
April 22. As a test of our code, we attempted to reproduce 
Drummond's associations by running our volu me search code 
on the same sample. Like lDrummon we often iden- 

tified overlapping associations in which one association was 
partially or wholly contained within another. We eliminated 
the smaller of any two associations that shared more than 66% 
of its members with the larger group. F ig. Q] shows that our 
result is slightly different from |D^mmondJ l20_QQj) although 
there is rough agreement in the distribution of the sizes of 
the associations. Table [2] compares the sizes and overlap be- 
tween the four larg est associations found in this work and in 
iDrummo'ndl ll2000l) . There is clearly a great deal of overlap be- 
tween the groupings found in both works but the relationship 
is not perfect. 

We believe there are three main reasons for the disagree- 
ment: (1) The orbital elements for the 708 NEOs are no 
longer identical to what they were five years ago. Most of 
these objects have much better orbits now due to new ob- 
servations provided to the MPC in the intervening time pe- 
riod; (2) The volume search technique is very sensitive to 
the input orbital elements. Convergence of the iterative tech- 
nique for identifying associations can alter membership by a 
few objects with changes in the Z?-criterion of only 0.01 (If 
all the difference between two orbits is in either a or e then 
Aa = 0.01 or Ae = 0.01 respectively). Reproducing the os- 
culating elements of the 708 NEOs used by Drummond would 
be a monumental task involving extracting only those detec- 
tions known to be associated with each object at the time that 
IDrummondl (120001) obtained the orbital element data and then 
determining the orbit; (3) The technique often identified over- 
lapping associations that shared many members. As described 
above, when the overlap between two groups exceeded 66% 
the smaller one was discarded which biases our method to- 
wards larger groups. 

Since we believe that our technique is demonstrably identi- 

2 We used Pikaia 1.2 [Charbonneau 2002), a FORTRAN version of a 
genetic search algorithm, as our non-linear least square fitting code to find 
the global minimum \ 2 - 
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cal to iDrummondl (120001) given the result shown in Fig.fflwe 
proceed with the use of the contemporary osculating elements 
for all 708 NEOs. 

3. STATISTICAL SIGNIFICANCE OF NEO GROUPS 

iDrummo'ndl (12000) tested his NEO grouping software on 
synthetic populations of NEOs but did not, apparently, take 
into account observational selection effects when creating 
those populations. This is not surprising considering the dif- 
ficulties involved in creating a realistic NEO population and 
then simulating the surveying capability of the many observa- 
tory programs that contributed to the 708 NEOs in his sample. 
Drummond tested the veracity of his groupings by performing 
the volume search on two sets of randomized orbital elements 
and comparing them to the results for the real set. He con- 
cluded that 50 to 70% of the groupings he had found might 
not be real. 

Since that time, a good four-dimensional orbital element 
and abs olute magnitude mode l of the NEOs has been devel- 
oped by iBottke et alJ ll2002al) and good matches have been 
obtained between simulated NEO surveys and the k nown ob- 
jects (e.g.. lRavmond et alJ2004llJedicke et al.l2003l) . 

To test the statistical significance of the NEO groups that 
we identified (see Fig. we generated 100 separate syn- 
th etic NEO models ac cording to the (a,e,i,H) distribution 
o flBottke et al.N2002al) . Each of the three angular elements of 
the orbits were generated randomly in the range [0, 2ir). We 
then passed e ach of them through a survey simulation simi- 
lar to those of lJedicke et alJ ( 120031) until the simulated survey 
detected 708 NEOs, identical to the size of the NEO sam- 
ple used in Drummond's (2000) analysis. Very briefly, the 
simulator attempts to mimic the average surveying capabil- 
ity of all surveys that have identified NEOs. It covers a 180° 
wide range in RA and from -30° to +80° in declination in the 
space of 14 days centered on new moon each lunation. We 
found that using a limiting magnitude of V = 20.5 (50% ef- 
ficiency) where the efficiency drops from 100% to 0% in the 
range from V = 19.75 to V = 21.25, and a minimum rate 
of motion to distinguish NEOs of 0.3°/day, yielded results 
that qualitatively agree with the observed distributions. The 
minimum and maximum range of the distributions of the 100 
synthetic populations is shown in Fig. [2] in which it is clear 
that our simulation nicely bounds the observed distribution of 
NEO orbital eleme nts and absolute magnitudes as utilized by 
IDrummondl (120001) . 

We used the volume search method to identify associations 
in our 100 synthetic NEO populations. In Fig. [3] we show 
the size distribution of all the associations detected in our 
synthetic NEO populations and the distribution of associa- 
ti ons found (by our c ode) in the same NEO sample as used 
in IDrummondl ll2000l) . Since no genetic families exist in the 
artificially generated NEO populations, the fact that we iden- 
tified associations containing more than 20 members in that 
data, and that the size distribution of the groupings for the 
real data is similar to that for the synthetic data sets, brings 
int o question the gen etic nature of the associations reported 
bv lDrummondl (120001) . 

Interestingly, both the synthetic and real data show a bi- 
modal distribution. When the size of the association is <15 
there is a power-law drop in the SFD. For associations con- 
taining >15 members the distribution is almost flat or mod- 
estly peaked near 22. We are unable to explain why the 
technique produces this bi-modal distribution. It is not the 
case that the synthetic model may contain synthetic genetic 



families because the actual NEO population from which it 
was derived contains genetic families. This is a result of 
the technique used t o generate the NEO model as detailed in 
IBottke et alJ ( 1200 2a). Briefly, the correlation between the or- 
bital elements of members in the actual NEO popula tion wa s 
lost in the fitting procedure employed by Bottke et all(l2002al) . 
When we generated our synthetic NEO population it therefore 
has no memory of any possible groupings in the real NEO 
population. 

We test whether the real NEO population is more clustered 
than the synthetic ones by determining if the two popula- 
tions are different or drawn from the same underlying dis- 
tribution function. The Kolmogorov-Smirnov (K-S) te st (e.g., 
von Mises 1 H%4l lKolmogorovlU933t iSmirnov lll933h is one 
of the most generally accepted techniques for testing the con- 
sistency of two one-dimensional distributions. It is applicable 
to non-gaussian distributions and works well on small sample 
sizes. We applied the K-S test to compare the size distribu- 
tion of the associations found in the real NEO population and 
that found in each of the 100 synthetic NEO populations. We 
found that the size distributions of the associations detected 
in 35 (54) out of the 100 synthetic populations are consistent 
with that of the associations detected in the real NEO popu- 
lation at a significance level 95% (85%). Figure [4] shows 
the cumulative distribution of the K-S probability that the syn- 
thetic data sets are consistent with the actual data. Thus, there 
appears to be litt le reason to believe that any of the associa- 
tions detected by IDrummondl (120001) in the actual NEO sam- 
ple are real. 

4. DECOHERENCE OF A SYNTHETIC A 1 -LIKE FAMILY 

In this section we determine the length of time that a syn- 
thetic genetic NEO family like Drummond's (2000) Al asso- 
ciation would be detectable using the Z?-criterion technique. 
This is important for understanding whether the Al associa- 
tion is real and also for calculating the expected number of 
NEO families that may be detectable in the future. 

A significant fraction of the Al association's mean orbit 
overlaps the MB so we assume that, if it is real, it was most 
likely created in a collision between a NEO and a MB aster- 
oid as opposed to one of the other four possible mechanisms 
identified in fJ2 Current computing techniques allow a re- 
alistic simulation of the entire process of the formation and 
evolution of aste roid families within a re asonabl e amou nt of 
clock time (e. e.. iMichel et alJl2002t fDurda et al.ll2004l) . We 
have perform ed a smooth-partic le hydrodynamics (SPH) sim- 
ulation (as in lDurda et all2004l) of a NEO-MB asteroid colli- 
sion from which an A 1 -like family was formed, tracked the re- 
accumulation process, and then followed the dynamical evo- 
lution of the individual f ragments. It should be noted that 
some recent works (e.g., Idell'Oro et al.ll2004l ICellino et alJ 
120041) suggest that detailed SPH models may underestimate 
the ejection velocities amongst the collision fragments. 

First, we specify the parent body for the collision: as- 
suming that the Al association is a genetic family we set 
the orbit of the synthetic Al parent body to the mean orbit 
of the actual 25 Al members, i.e., (a,e,i,Q,uj) = (2.21 
AU, 0.49, 4.1°, 194.3°, 123.3°). 

The lower limit on the size of the parent body is the sum 
of the volume of the 25 A 1 me mbers which i s dom inated by 
the largest fragment, 1627 Ivar. IVeeder et al.l (119891) reported 
radiometric observations of this obj ect and 2368 Beltr ovata 
(the third largest member of Al) and lDelbo et alJ ( 120031) pro- 
vide independent radiometric observations of 1627 Ivar. De- 
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pending on the thermal model used, IVeeder et alJ (119891) de- 
termined a visual geometric albedo of 0. 08 or 0.12 t o 1627 
Ivar and 0.05 or 0.13 to 2368 Beltrovata. iDelbo et all d2003l) 
found a visual geometric albedo for 1627 Ivar of 0.20, 0.15, 
or 0.050 depending on whether the STM, NEATM, or FRM 
thermal model was used. On the other hand, the effective di- 
ameter for 162 7 Ivar was estimate d to be 8.5±3 km from radar 
measurements dOstro et al.ll 19901) . Given the widely varying 
albedo estimates from radiometric observations and thermal 
modelling, we decided to use the radar measurements as the 
basis for our constraints on, and best estimate for, the diame- 
ter of the Al parent body. 

Table[2shows the absolute magnitude {H), taxonomic type 
(where available), and effective diameter (lower bound, best 
estimate, upper bound) for each object in Al. The H val- 
ues for each object are from the Lowell Observatory's AS- 
TORB database (ftp://ftp.lowell.edu/pub/elgb/astorb.htmli. 
The lower limit for the Al parent body is shown for three 
different albedo values at the bottom of the table. Note that 
the albedo values (pv = 0.07 — > 0.306) in this table encom- 
pass all of the albedo estimates from radiometric measure- 
ments (0.07 — > 0.26) for 1627 Ivar, hence we believe 9.8±3.5 
km is a suitable estimate for the size of the Al precursor. For 
the purpose of our simulations we round this number to a syn- 
thetic Al prec ursor diameter of 10 km. 

We note that Tanga etal] (119991) claim that a better prelim- 
inary estimate of the parent body size based on geometrical 
considerations is obtained merely by summing the diameters 
of the first and third largest asteroids in the family. They state 
that this is a much better estimate of the parent body diam- 
eter when the volume method (as used in this analysis and 
described above) yields a diameter less than the sum of the 
two largest family members. Using their technique our best 
estimate for the diameter of the Al association parent body 
would be about 12 km. 

With the target asteroid specified we then determine the rel- 
evant properties of the impactor. We assumed that the orbital 
elements of the impactor are equal to those of a MB object that 
is likely to collide with the Al precursor. We identified these 
orbital elements in the following manner. The MB asteroids 
are almost complete to H <~ 14.5, so we made a histogram 
of the semi-major axis of all MB objects with H < 14.5 that 
could hit the Al precursor. Since the Al precursor's heliocen- 
tric distance ranges from perihelion at ~ 1.1 AU to aphelion 
at ~ 3.3 AU, we examined all MB semi-major axes with 1.1 
AU < a(l — e) < 3.3 AU and found that the most prob- 
able semi-major axis for the projectile is ~ 3.1 AU. Simi- 
larly, the most probable eccentricity is ~ 0. 1 1 and inclination 
is ~ 10.0°. Therefore, we select as our impactor an object 
with (a, e, i) = (3.1 AU, 0.11, 10.0°). 

We assume that the collision takes place at a speed equal 
to the most probable collision speed for an object with Al's 
mean orbit — 8.36 km/s (Bottke 2003, personal communica- 
tion). 

The most likely location of an impact event is near the aphe- 
lion of the Al precursor since it spends most of its time near 
that heliocentric distance. The range of overlapping heliocen- 
tric distances between the orbits of the Al precursor and the 
projectile is from 2.94 AU to 3.30 AU, but the actual range 
where the impact could happen is more restricted because we 
have fixed the impact speed. At higher heliocentric distances 
(> 3.09 AU) the relative impact speed dictates that the relative 
inclination (the angle between the planes of the two orbits) 
must be larger than the highest possible inclination difference, 



which is the sum of the two inclinations (14.1°). The loop in 
Fig.|5]shows the relationship between the relative inclination 
and the heliocentric distance of the impact when constrained 
by our choice of the impact speed. The dashed lines in the 
figure bound the region of permitted relative inclinations be- 
tween the Al precursor and the projectile, i.e., smaller than 
i\ + 22 and larger than | i\ — | where the subscripts 1 and 
2 refer to the target and projectile. For our simulations we 
used the smallest and largest allowable heliocentric distances 
from Fig. 0(2.94 AU and 3.09 AU) and the one at which the 
relative inclination is the smallest (2.98 AU) in order to inves- 
tigate the dependence of the decoherence time scale upon the 
impact location. 

Since we require that the collision occurs at a specific helio- 
centric distance and speed, the three angular orbital elements 
of the impactor can be determined because the location of the 
Al precursor at the time of collision is known and the im- 
pactor has no choice but to cross the same point at the same 
time. We note that there are two solutions — the impactor 
may collide on either hemisphere relative to the velocity vec- 
tors of the two objects. We present the results for only one 
case since both give essentially identical results. 

We modelled the initial stag es of the impact with the 3- 
dimensional SPH code SPH3D dBenz and Asphaug 119951) . In 
keeping with our pattern of favoring the most probable event, 
we selected the smallest possible impacting object that would 
produce a barely catastrophic disruption. Based on our previ- 
ous experience with impact simulations, and scaling the val- 
ues for the velocity of the impact and size of the target, we 
selected an impactor of 0.25 km diameter. We will argue later 
that this is probably too large an impactor but we believe that 
the results of the simulation are still of practical utility. Sim- 
ilarly, the most probable impact angle is 45° as used in this 
simulation. 

The 10 km diameter target was comprised of 100,000 SPH 
particles and, to roughly match the SPH particle volume den- 
sity of the target, the impactor is comprised of 700 particles. 
In this simulation the Al precursor was catastrophically dis- 
rupted and shattered into thousands of independent fragments 
by the impactor. The SPH phase of the simulation followed 
the impact for 100 seconds after which time the ejecta flow 
fields were well established and no further fragmentation or 
damage was observed. 

The outcome of the SPH simulation was handed off as the 
initial conditions for an N-body si mulation using the cos- 
mologically derived pkdgrav code dRichardson et al \ [ 20001 
ILeinhardt etail 120001: ILeinhardt and Richardsonl 120021) that 
followed the trajectories of the ejecta fragments for sufficient 
time to allow gravitational re-accumulation of th e larger col- 
lision fragments. This technique was used bv iDurda et all 
d200l to model the formation of asteroid satellites resulting 
from large impacts on 100-km scale MB asteroids. With the 
small sized target assumed here for the Al precursor, there 
is not much gravitational re-accumulation among the largest 
fragments and the N-body phase of the simulation need only 
be run for 200 N-body timesteps (each step corresponding to 
about 50 s) in order to establish fragment positions and ejec- 
tion velocities. 

From the velocity vectors of the Al precursor and the im- 
pactor at the collision location we can derive the three Eu- 
ler angles between the impact space, where the simulation 
was performed, and the external heliocentric orbital space in 
which we follow the post-collision dynamical evolution of 
the fragments. The initial orbits of the re-accumulated frag- 
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ments are handed over to Mercury6 (IChambers 119991) to fol- 
low their dynamical evolution. The SPH simulation produced 
1243 fragments down to a diameter of 216 m (pre-determined 
by the resolution of the simulation). Note that the smallest 
diameter amongst known members of the real Al association 
is only 130 ± 50 m (see TableQ. 

Only a fraction of the fragments would be discovered in re- 
ality, so we applied an ad hoc observational selection function 
to determine which of the synthetic fragments were detected. 
The correction factor is a complicated function of a survey's 
performance characteristics and the orbital ele ment distribu- 
tion of the detected objects iTedicke et al.l2002l) . Since all the 
synthetic Al family members have essentially identical or- 
bital elements we would like to determine the selection effect 
as a function of H for all the surveys that contributed to find- 
ing NEOs with orbital elements similar to the Al mean orbit. 
This would be difficult if not impossible to determine in prac- 
tice. Instead, we chose simply to divide the set of all known 
708 NEOs at the time of Drummond's (2000) analysis as a 
function of absolute magnitude, n( H), by the debiased " true" 
distribution of NEOs according to Bottke et al. (2002a|). We 
fit the resulting distribution to a function of the form: 

P(H) = -ff=J , (8) 

e -tr- + 1 

and found L = 15.82 ± 0.19, W = 1.15 ± 0.11. P(H) 
gives the probability that an NEO with absolute magnitude 
H would have been discovered by the time of Drummond's 
(2000) analysis. 

We applied this discovery probability with the central val- 
ues for L and W to the re-accumulated synthetic Al family. 
An albedo of 0.128 (the same as was used for the best esti- 
mate of the size of the Al precursor as in Tabled was used 
to convert the size of the objects into an absolute magnitude. 
Running many trials for this process (which objects "survive" 
the selection function of equation[8]is random) yielded a mean 
size of an observed synthetic Al family of 25 ± 4 members 
compare d to the 25 member s of the actual Al association re- 
ported bv lDrummon 51 II2000I The size distribution of our syn- 
theticpopulation is compared to the actual size distribution in 
Fig. E| 

It appears that our model does not reproduce the actual Al 
association's size distribution. This could be due to our hav- 
ing selected an impactor that caused far too much damage to 
the target, shattering it with such violence that the fragments 
were unable to gravitationally accumulate into larger objects. 
However, considering that our dynamical integrations of the 
fragments consider only gravitational effects, the SFD of the 
fragments is probably not important to the remainder of this 
study. 

We now examine the decoherence of the synthetic Al fam- 
ily as its members dynamically evolve under the gravitational 
influence of the Sun and eight planets (Mercury through Nep- 
tune). It was only necessary to integrate the motion of the 25 
fragments that were detected rather than the ensemble of 1242 
fragments. The orbits were propagated forward in time us- 
ing a hybrid symplectic/Bulirsch-Stoer integator (Mercury6, 
Chambers 1999) with a step-size, r, of 30 days (i.e., 1/40 
of the shortest initial orbital period amongst the fragments). 
We verified the result by comparing it to an identical run on 
the same fragments using the general Bulirsch-Stoer integra- 
tor with the same step-size. 

In reality, this evolution occurs within a population of back- 
ground objects that are also interacting gravitationally with 



the major objects in the solar system and with each other. In 
order to save time, for the purpose of this study, we decided 
not to integrate the motion of all the background objects. In- 
stead, we used two different static background NEO popula- 
tions: Bl was a randomly selected synthetic NEO population 
(from with 708 objects and B2 was formed using the same 
708 NEOs as used in Drummond's study but removing the 25 
actual Al members. This means that our static background 
populations possessed a slightly different mean density in or- 
bital element space with a total of 733=708+25 objects in the 
Bl case and only 708 objects in the B2 case. This introduces 
a difference in the mean density of only a few percent. It 
would be difficult to correct because it would require selec- 
tively removing 25 members from the synthetic Bl popula- 
tion that somehow match the original Al members. Instead, 
we chose to ignore this small difference in density and note 
that, in any event, it makes the identification of groupings in 
the synthetic B 1 data slightly more likely. We introduced our 
synthetic Al family members into the background NEO pop- 
ulations and search for groupings every 2,000 years as the Al 
family evolves in time. Since our purpose was to track the 
evolution of the Al family, we only performed the volume 
search for candidate groupings around the 25 detected frag- 
ments as opposed to searching the entire sample. The mean 
orbit of the resulting association may be beyond the boundary 
of the original fragments because interloping asteroids may 
be incorporated into the final association due to the iterative 
aspect of the volume search technique. 

While our d efinition of associa tions, strings and pairs 
matches that of iDrummon our manner of selecting 

the final groups may be different. We searched for associa- 
tions using the (known) synthetic Al members as seeds and 
select the association that includes the most members of the 
original synthetic family. If multiple associations meet this 
criteria then we selected the association containing the small- 
est number of interlopers. Similarly, strings were identified 
using the synthetic Al members as seeds and selecting the 
string that contains the most members of the original associ- 
ation. Once again, if multiple strings met this criterion we 
used only the string that contained the least interlopers. Fi- 
nally, pairs were identified only for objects within the original 
synthetic Al association in the current timestep. 

Fig- demonstrates the decoherence of an Al-like family 
from creation until it is 1 Myr old. The results for other 
choices of the 25 detected fragments, the NEO background 
and for the heliocentric distance of impact appear nearly iden- 
tical. It shows that the rapid orbit evolution of NEOs mixes 
the genetically related synthetic Al objects and background 
objects within ^300 Kyr so that the association becomes un- 
detectable to this technique. The fact that all our runs have 
similar results indicates that the decoherence time scale is in- 
dependent of the choice of the background NEO population 
and the location where the precursor is disrupted. 

5. DECOHERENCE OF THE ACTUAL Al ASSOCIATION 

TheAl associat ion is the largest among the 14 identified by 
iDrummondl (120001) and all its 25 members are included in the 
largest association identified in this study (see Table|2j. It is 
the most statistically significant association found using our 
technique and the 25 common members have a mean orbit of 
(a,e,i) ~ (2.21 AU, 0.49, 4.1°). In g]we argued that the 
Al association must have derived from a parent body > 10 
km in diameter if it represents a genetic NEO family, yet we 
calculated in Sj^that the production of a family from the dis- 
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ruption of an object this large is extremely unlikely. However, 
just because an event is unlikely doesn't mean that it did not 
occur. 

In this section we use two different methods to examine the 
possibility that Drummond's (2000) Al association is a bona 
fide genetic family or due merely to chance alignments or ob- 
servational biases. First, we integrate forward the orbits of all 
the Al members and see how long it remains detectable as an 
association. If the association is real we would expect that the 
members remain an association for some time into the future. 
Second, we integrate the orbits backward in time to see if their 
angular orbital elements converge, indicating a common point 
of origin at some time in the past. The backward integration 
is justified for this association because we know that, if it is 
a genetic family, it can not be very old. The latter technique 
has been used bv lNesvornv et alJ ll2002l) to determine the age 
of the Karin cluster within the MB. 

The orbits of the 2 5 me mbers in Al were propagated for- 
ward (and backward, i )5.1> in time using the general Bulirsch- 
Stoer integrator with r = 20 days (i.e., 1/40 of the short- 
est initial orbit period amongst the actual Al members). 
The current orbital elements of the Al members are ex- 
tracted from the MPCORB database as of 2005 March 1 
(ftp://cfa-ftp.harvard.edu/pub/MPCORB/! that contains accu- 
rate orbital elements derived from multi-opposition observa- 
tions of the NEOs. Once again, the integration took into ac- 
count the gravitational perturbations due to all major plan- 
ets (Mercury through Neptune) but neglected the evolution of 
background NEOs. Every 2,000 years we searched for group- 
in gs (pairs, strings & a ssociations) in the same sample as used 
by iDrummondl ( 120001) using the techniques described above 
but replacing the 25 Al members with their evolved orbits. 

For the same reasons as in we on ly searched for den- 
sity enhancements above the local background using the 25 
Al members as seeds in the volume search. We searched for 
strings and pairs within the Al association at each time step 
in a manner identical to that described in f@] 

Fig.|8]shows the evolution of Al during the next 500 Kyr. 
The Al members rapidly evolve away from each other as 
shown in Fig.|HJ\ — the average D-criterion amongst all Al 
members increases from 0.2 to >0.6 in ~20 Kyr and to nearly 
1.0 after slightly more than 200 Kyr. At the first time step 
{i.e., now), all the 25 Al members are included in our 29- 
member association (Fig. HJi), but after just two time steps 
(4,000 years) the association is essentially undetectable. In 
the following 500 Kyr the largest A 1 -like associations contain 
anywhere from a few to about 30 objects but those groups con- 
tain only a few of the original Al members. In other words, 
those associations are quickly and heavily contaminated by 
interlopers. 

The situation is similarly bleak for strings (Fig. |8jZ) and 
pairs (Fig. |8jD) detected within the associations at each time 
step. The largest string detected at the first timestep contains 
only 5 members and after just 20,000 years not a single time 
step ever contains more than 4 members of the original Al 
cast. Similar behavior is observed for the pair counts within 
the original Al members. With 25 objects there are 300 possi- 
ble unique pairings but the maximum number of pairs within 
the Al association is 10 for only the very first time step. 

The result that the Al association is essentially unde- 
tectable within ~10 Kyr suggests that it too is merely a sta- 
tistical fluctuation in the density of NEO orbital elements and 
argues against the genetic nature of Al. It is extremely un- 
likely that such a large association would be detected just as 



this technique lost its ability to identify it due to the orbital 
decoherence of its members. 

5.1. Evolution oftheAl association's Q and lo 

iNesvornv et alJ ll2002ft integrated the orbits of members of 
the Karin cluster backwards in time to identify when the lon- 
gitude of ascending node (17) and argument of perihelion (cu) 
of all members were tightly clustered. These angular orbital 
elements evolve rapidly and secularly under the influence of 
the gravitational perturbations of the planets. Determining the 
point in time at which the values were tightly clustered al- 
lowed them to accurately date the age of that small and new 
asteroid family. The evolution of these angular elements will 
be considerably faster for NEOs but we tested for their con- 
vergence despite the low probability of success. 

Figure|9]shows the evolution of the node and perihelion an- 
gles for the last 50 Kyr. We only show the evolution over this 
reduced time period because it is representative of the entire 
evolutionary period of this association (see next paragraph for 
discussion). Note that at t = the values are not consistent 
with being randomly distributed since the D-criterion tech- 
nique specifically selects sets of objects with non-random or- 
bits (in £1 and lu as well as a, e and i). All 25 values of the 
longitude of ascending node are distributed in the restricted 
range 125° < f2 < 325° and while the longitude of perihe- 
lion is similarly restricted in the range 0° < fl < 180°. 

To evaluate whether these orbital elements become less ran- 
dom in the past, in other words, more "clumped", we deter- 
mined the probability that the elements were consistent with 
being flatly distributed using the K-S test at each time step. 
This test was uninformative, the distribution changing rapidly 
and randomly from consistent to inconsistent over the entire 
300 Kyr period. The problem is that there are many ways for 
a distribution to be "clumped". Instead, we resorted to exam- 
ining the distribution of £1 and u in terms of their range, and 
RMS spread within that range, as a function of time as shown 
in Fig. \lQ The range of the two elements and their RMS 
spread within that range is at a minimum at the present time. 
In other words, at no time in the past were these angular ele- 
ments more "clumped" than at the present time. The "dumpi- 
ness" at the present time is easily explained as a consequence 
of the D-criterion technique selecting clumps in their orbit 
distribution. Once again, this result contradicts the hypothe- 
sis that the Al association represents a genetically linked set 
of asteroids. 

6. IDENTIFYING ACTUAL GENETIC FAMILIES 

The fact that the volume search method identified a num- 
ber of "associations" in synthetic NEO populations where no 
genetic families exist (||3j indicates that the method, by it- 
self, can not distinguish between genetic families and random, 
over-dense regions in the orbit element space. In this section 
we show that, in principle, it should be possible to separate 
the real genetic families from the background by combining 
results from the volume search and string/pair searches. 

In we showed that a genetic NEO family loses D- 
criterion coherence ~ 250 Kyr after creation. However, 
even after the family diffuses into the background the vol- 
ume search method still detects associations with > 15 ob- 
jects (Fig.Qi) but most of the members of these associations 
are interlopers from the background. The question is then how 
to distinguish between associations containing actual geneti- 
cally related members and those containing interlopers. 
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Comparing Fig. 03 with Figs. reveals a clear dif- 

ference between the evolution of the association size and the 
string size or total number of identified pairs — after the fam- 
ily loses coherence the association size becomes very noisy, 
jumping erratically between small and large families; while 
the string size and number of pairs drops to a near steady- 
state value with little variation from time step to time step. 
We utilize the difference in behaviour between associations, 
strings and pairs to differentiate between real and false ge- 
netic NEO families in Fig. ^2 in this simulation the strings 
are found using the members of the association identified at 
each time step as seeds. The pairs are identified only amongst 
the members of the association at each time step and the av- 
erage D-criterion is also calculated only amongst the same 
members. Our intent is to mimic the fact that when searching 
for new genetic NEO families the members of the association 
will not be known a priori as in Figs.Qand|8] In that figure 
the pair fraction is the total number of pairs identified within 
the association divided by = n(n — l)/2, the total number 
of possible unique pairs within the association when n is the 
number of objects it contains. The average D-criterion value 
is the sum of the D-criterion values between each unique pair 
of objects in the association also divided by C\. These two 
parameters essentially measure the "tightness" of an associa- 
tion. We expect that fresh genetic NEO families will be very 
"tight" with pair fractions approaching unity and average D- 
criterion value approaching zero. There is a very strong corre- 
lation between the pair fraction and average D-criterion value. 

The data in Fig. ^2 are f° r the synthetic Al family of f|4] 
Each data point represents a single time step in the dynamical 
evolution of the synthetic Al family with the different sym- 
bols representing time since its formation. There is clearly 
a strong time-dependent evolution that may be used to dif- 
ferentiate between real and false genetic NEO families. The 
0.1 D cuto f f is found epirically to yield the clearest separation 
between associations detected at different epochs. The black 
crosses in the figure represent the results of the same analy- 
sis applied to our synthetic NEO background populations (see 
The groupings found in those populations (with no fami- 
lies) overlap with the ones found after the decoherence of the 
synthetic Al family but are well-separated from the groupings 
found before the decoherence time of 250 Kyr. The current lo- 
cation of Drummond's (2000) Al association is indicated in 
Fig. II llwith a red star, and the other associations (A2 - A14) 
are marked with blue triangles. 

7. DISCUSSION 

We believe that the results described above bring into se- 
rious que stion the reality of any of the NEO associations re- 
ported bv lDrummon l m6dt) . It is unlikely that a large NEO, 
similar in size to the precursor to Drummond's (2000) Al as- 
sociation, would have disrupted during the last 300 Kyr, the 
time frame during which the family of asteroids would be de- 
tectable with the D-criterion orbital similarity technique. The 
fact that th e largest, and presu mably most likely, association 
reported b v lDrummondl (12000 ) would be unidentifiable in just 
a few hundred thousand years implies that, if it is real, it was 
found just in the ~250 Kyr period before the association be- 
comes totally decoherent (see Figure 0. This represents a 
mere ~10% of the dynamical lifetime of typical NEOs and 
is in agreement with the decoh erence time scales for NEO 
meteorite streams as reported bv lPauls and GladmanN2005l) . 
Furthermore, the SFD of associations detected in the actual 
NEO sample matches the same distribution for associations 



detected in a random sample of NEOs that accounts for ob- 
servational selection effects. 

Confusing the issue slightly, we applied our technique to 
all of Drummond's (2000) original associations but using his 
more relaxed (D cuto f f) threshold on the D-criterion of 0.1 15 
instead of our value of 0. 1 . All but one of these associations 
remained firmly in the region of random background fluctua- 
tions shown in Fig. ^2 Only when using the original members 
of the A14 association as seeds with the relaxed Z?-criterion 
threshold did we detect a 5-member string (as opposed to a 
2-member "string" with D cuto ff = 0.1) containing all four 
original A14 members (String Size/Association Size = 1.2). 
At the same time, we found 3 pairs amongst the four A14 
members (i.e., Pair Fraction = 0.5) compared to only 1 pair 
when D cuto ff =0.1. In this case, A14 would be located 
in the more interesting region of Fig. where genetic fam- 
ilies may be identified, well beyond the noise region defined 
by the fake detections within the synthetic NEO populations. 
This hints that A14 may be a genetic NEO family but since 
it has only 4 members, losing or gaining a single member in 
the string would cause a 25% difference in the ratio of the 
string size to association size. Furthermore, our diagnostic 
technique is based on the simulation of a 25 member A 1 -like 
family rather than a 4 member A14-like family so the com- 
parison is not strictly applicable. 

On the other hand, we show in Section|4]that a NEO fam- 
ily can maintain coherence in it's orbital elements for >200 
Kyr. So the Z?-criterion technique may be useful for identify- 
ing NEO families that must then undergo further study before 
being classified as real. The supplementary studies should 
include at least a statistical study of whether groups of the 
identified size are unlikely in a similarly sized synthetic data 
set that incorporates all the known observational selection ef- 
fects. 

Strong support for the veracity of a putative association 
found by the volume search method may be obtained from 
its location in diagnostic diagrams similar to Fig. ^3 All the 
points with pair fraction > 0.2 and string size / association 
size > 0.8 correspond to positive detections of real genetic 
associations with very few interlopers. Alternatively, the fig- 
ure may be considered as demonstrating the intrinsic limits of 
Z?-criterion based methods in identifying families - only the 
tightest families (i.e., youngest ones) can be solidly detected. 
We note that the current location of Drummond's (2000) Al 
is solidly in the territory occupied by random fluctuations. 

We understand that we h ave only examined in detail a sin- 
gle association identified bv lDrummondl J2000I) . It is possible 
that parent-body NEOs in other orbits may be more likely to 
disrupt and form longer living, more easily identified NEO 
families. It is also possible that one of the other methods of 
producing an NEO family may produce longer lived or more 
easily identifiable families. The different production mecha- 
nisms may even show a difference in the distribution of the 
orbital elements of their members during the early stages of 
their separated evolution. 

For instance, the tidal disruption process (e.g., 
iRichardson et al.lll998l) seems capable of producing contact 
and close binary asteroids and also families of objects on ver y 
similar orbits (e.g., Comet Shoemaker-Levy 9. lBossl dl994l) : 
crater chain s on t he various planets and their satellites, 
iBottke et alJ J1997I) - ). The fragments from this production 
process may well be detectable as a family long after the 
planetary close approach that produces the chain of objects. 
We plan on studying this possibility in future simulations. 
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The dynamical simulations of both the actual and synthetic 
Al families (i ]5]and $\ did n ot account for the Yarkovsky ef- 
fect (see iBottke et alTl2002bl for a review) that is known to 
cause diffusion of orbital elements for small objects in the 
MB. The effect of this non-dynamical force will be corre- 
spondingly greater for NEOs due to their approaching much 
closer to the Sun. Even so, the expected drift rate for a 100 
m diameter NEO with a = 1 AU, i = 45° and typical den- 
sity, albedo and surface conductivity, is only ~2x 10~ 4 AU 
in 10 5 years (Nesvorny, personal communication). Thus, the 
Yarkovsky effect operates far too slowly to contribute to the 
decoherence of NEO families. 

It is clear that contemporary studies of genetic NEO fami- 
lies are hampered by the small number of known NEOs thus, 
larger and deeper surveys are needed. Pan-STARRS will be 
the next major survey to come online and will detect as many 
asteroids and comets in one lunation as are currently known. 
Pan-STARRS will revolutionize our understanding of NEOs 
— while 3,319 NEOs are known as of April 2004 there are 
hundreds of thousands of NEOs brighter than Pan-STARRS' 
expected limit of R ~ 24 mag. 10,000 new NEOs should be 
discovered in the first year of operations alone. By the end of 
the first year of Pan-STARRS operation, the completeness of 
NEOs of 1 km (100 m) in diameter will be boosted to at least 
90% (20%) from the current 75% (5%), and after 10 years of 
operation Pan-STARRS will discover almost all of the NEOs 
larger than 1 km diameter and more than half of the NEOs 
bigger than 100 m diameter. With a good understanding of 
the observational selection effects in a single survey it should 
be possible to identify NEO families or determine that there 
are none and use this information to constrain the collision 
rate of NEOs. 

8. CONCLUSION 

We find it unlikely that any of the enhancements in orbital 
element space for the known NEO population are due to a ge- 
netic relationship between the member objects. The primary 



motivation for this conclusion is the fact that synthetic data 
sets that incorporate realistic distributions of NEOs and obser- 
vational selection effects show the same distribution of den- 
sity enhancements as the actual population. We also base our 
conclusion on the fact that the current understanding of NEO 
collision rates makes it extremely unlikely that a large enough 
NEO could have disrupted recently enough to allow its frag- 
ments to maintain enough coherence to remain detectable by 
the D-criterion technique. 

We have identified a new technique that will allow future 
searches for genetic families within the known NEO popu- 
lation. The method relies on the tight clustering in the D- 
criterion during the first couple hundred thousand years after 
the production of the family. Once an association is identified, 
it is searched for strings and pairs. If the number of found 
pairs is more than ~20% of the maximum possible number of 
pairs in the association and the ratio of the maximum string 
size to association size is > 0.8, then the association is likely 
to be a real genetic family. 

Upcoming improvements in NEO survey technology such 
as Pan-STARRS may provide a large enough sample of NEOs 
in a single, well-characterized survey, to allow identification 
of NEO families or set a limit on their number. This, in turn, 
will allow dynamicists to refine their models of the collisional 
evolution of the solar system. 
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TABLE 1 

The A 1 association of IDrummondI j2000D in order of increasing absolute magnitude (H). This is the largest 

ASSOCIATION IDENTIFIED WITHIN THE NEO POPULATION KNOWN AT THE TIME. OUR NOMINAL VALUES FOR THE ASSUMED 
DIAMETERS OF THE OBJECT CORRESPOND TO THE INTERMEDIATE ALBEDO VALUE OF p v = 0.128. 
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TABLE 2 

A COMPARISON BETWEEN THE FOUR LARGEST ASSOCIATIONS IDENTIFIED IN THIS ANALYSIS AND In|DRUMMOND I 200( 3). TH E FOUR LARGEST 
GROUPS FOUND IN THIS ANALYSIS CONTAIN FOUR OF THE ASSOCIATIONS IDENTIFIED BYrDRUMMONp UiOOOIl . 
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FIG. 1 . — Size distribution of associations published in Drummond 1 2000) (solid line) and the associations detected by the volume search code of this analysis 
(dashed line). 



Identifying NEO Families 



13 





Inclination (deg) Absolute Magnitude 



FIG. 2. — Minimum and maximum (lower and upper dashed lines respectively) distribution of 100 synthetic observed NEO populations designed to model the 
distribution of NEOs known at the time of Drummond's (2000) study (solid line). A) semi-major axis B) eccentricity C) inclination D) absolute magnitude. 
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Association Size 



FIG. 3. — Average size distribution of associations detected from 100 synthetic NEO populations (solid line) and the distribution of associations found in this 
analysis in the same sample as used in Drummon 3 HxIJ (dashed line). 




FIG. 4. — Cumulative distribution of the Kolmogorov-Smirnov probability that the size distribution of associations found in a synthetic NEO population and 
the one from the actual NEO population are consistent. For example, at a confidence level >90%, the size distributions of associations detected in ~40 of the 
100 synthetic NEO populations are statistically consistent with the same distribution found in the actual NEO population. 
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Heliocentric Impact Distance (AU) 

FIG. 5. — Relative inclination vs. heliocentric impact distance for the Al precursor and projectile (see text for details). The loop shows the relative inclinations 
required in order to provide the relative impact speed (8.36 km/s) as a function of the heliocentric impact distance. The area bracketed by the dashed lines 
indicates the region of possible relative inclinations between the Al precursor and the projectile, i.e., smaller than the sum of and bigger than the difference 
between their orbit inclinations. 




FIG. 6. — Cumulative size distribution of the actual Al association (dots) and the synthetic Al-like family before (thin line) and after (thick line) applying an 
ad hoc observational selection effect. 
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FIG . 7 . — Evolution of a synthetic A 1 -like association in the next 1 Myr. These results use the B 1 background population and an impact heliocentric distance of 
2.98 AU for the collision (see text for details). A) the mean D-criterion value among all synthetic Al members, B) the largest detectable Al-like association (thin 



tring 



0.1 



line) and the number of objects that are original members of the synthetic family (thick line) in each detected group, C) the largest string with D, 
(thin line) and the number of objects that are original members of Al (thick line), D) total number of all pairs within the original synthetic Al members meeting 



the D„ 



0.1 threshold. 
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FIG. 8. — Evolution of Drummond's (2000) Al association in the next 500 Kyr — The sub-figures are otherwise identical in description to those in Fig.lTlbut 
are reproduced here for clarity. A) the mean D-criterion value among all original Al members, B) the largest detectable Al-like association (thin line) and the 
number of objects that are original Al members detected in each detected group (thick line), C) the largest string with D string = 0.1 (thin line) and the number 
of objects that are members of the original Al association (thick line), D) total number of all pairs amongst the original Al members meeting the D pa i r =0.1 
threshold. There are a total of 300 (Cf 5 = 25 X 24/2) possible pairs within the 25-member Al association. 
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FIG. 9. — Time evolution of the longitude of the ascending node, Q (top), and of the argument of perihelion, u> (bottom), for the 25 members of Drummond's 
(2000) Al association. The orbital elements are shown relative to the value for 1998 QQ52 (that appears as the straight line at 0°) because its orbit is the closest 
to Al's mean orbit. 
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FIG. 10. — Left - evolution of the range (max - min) of f2 (dotted) and w (solid) in the last 500 Kyr for Drummond's (2000) Al association. Right - evolution 
of the RMS spread of the same angular elements. In creating these figures we have taken into account the wrap-around at the 360° — > boundary. 
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FIG. 1 1 . — Left - pair fraction vs. string and association size ratio diagram. Right - average Z)-criterion value vs. string and association size ratio diagram. 
Large magenta solid dots — associations detected between T = and 210 Kyr; large magenta circles — associations detected between T = 210 Kyr and 250 Kyr; 
small magenta dots — associations detected after 250 Kyr; black crosses - associations detected in the synthetic NEO populations; red star — current location of 
Drummond's (2000) A l association; blue open triangles — current locations of Drummond's (2000) associations A2 - A14. 



